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ABSTRACT 

Gravity-capillary  standing  waves  in  water  of  arbitrary  uniform  depth  are 
considered.  The  classical  perturbation  calculation  yields  unbounded 
coefficients  for  some  critical  values  of  the  depth.  A  perturbation  solution 
valid  at  the  first  critical  value  of  the  depth  is  derived.  It  is  found  that 
two  solutions  exist  at  this  critical  value.  Numerical  computations  indicate 
that  these  solutions  are  members  of  two  different  families  of  solutions. 
Graphs  of  the  results  are  Included. 
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SIGNIFICANCE  AND  EXPLANATION 


In  recent  years  Important  progress  ham  b«an  achieved  in  tha  undar standing 
of  tha  affaot  of  surfaca  tansion  on  nonlinmar  fraa  surfaca  flow  pr obi  mas.  For 
axampla  Schwarts  and  Vanden-Broeck  (1979)  constructed  solutions  of  high 
accuracy  for  gravity-capillary  progressive  waves.  Their  results  indicate  the 
existence  of  a  number  of  different  continuous  families  of  solutions. 

In  the  present  paper  we  consider  gravity-capillary  standing  waves  in 
water  of  arbitrary  uniform  depth.  This  problem  was  first  considered  by  Concua 
(1962).  He  calculated  the  solution  to  third  order  as  a  power  series  expansion 
in  the  wave  amplitude.  He  found  that  some  of  the  series  coefficients  are 
unbounded  for  some  critical  values  of  the  depth. 

We  present  a  perturbation  solution  valid  at  the  first  critical  value  of 
the  depth.  We  show  that  two  solutions  exist  at  this  critical  value.  In 
addition  we  use  the  numerical  scheme  derived  by  Vanden-Broeck  and  Schwarts 
(1981)  to  compute  the  solution  in  the  neighborhood  of  the  first  critical  value 
of  the  depth.  We  show  that  the  two  solutions  obtained  at  the  critical  value 
are  members  of  two  different  families  of  solutions.  Similar  properties  were 
found  by  Schwartz  and  Vanden-Broeck  (1979)  for  gravity-capillary  progressive 
waves  in  the  neighborhood  of  the  first  critical  value  of  the  capillary  number. 


The  responsibility  for  the  wording  and  views  expressed  in  this  descriptive 
nmsiry  lies  with  MXC,  and  not  with  the  author  of  this  report. 


MONLXNBAR  GRAVITY-CAPILIARY  STANDING  NAVSS 
XN  MATIN  ON  ANBXTNANT  UMXFOMN  DCTTH 


Jaan-Marc  Vandan-Broeck 


1.  INTRODUCTION 

The  problem  of  gravity  standing  waves  in  watar  of  arbitrary  uniform  depth 
waa  solved  to  third  order  by  Tadjbakhsh  and  Keller  (1960).  Their  method  was 
applied  by  Concus  (1962)  to  solve  the  sere  general  problem  which  includes 
capillary  as  well  as  gravitational  forces. 

These  perturbation  expansions  were  obtained  by  imposing  a  uniqueness 
condition  which  excludes  certain  fluid  depths.  Concus  (1964)  showed  that  the 
values  of  the  depth  excluded  by  this  condition  fora  a  denuaerably  infinite  set 
which  is  densely  distributed  over  the  entire  positive  real  line.  Xt  is 
therefore  essentially  iapoasible  to  satisfy  the  uniqueness  condition  in 
practice.  However  the  solution  obtained  by  Tadjbakhsh  and  Xellar  (1960)  is 
satisfactory  since  it  is  defined  for  any  value  of  the  depth  including  those 
excluded  by  the  uniqueness  condition.  These  results  are  confirmed  by  the 
numerical  calculations  of  Vanden-Broeck  and  Schwarts  (1981). 

The  use  of  the  uniqueness  condition  in  the  general  problem  with  surface 
tension  results  in  unbounded  series  coefficients  for  certain  values  of  the 
depth  (Concus  (1962)).  Although  these  values  of  the  depth  were  excluded  by 
the  uniqueness  condition,  the  perturbation  solution  is  clearly  not 
satisfactory  for  values  of  the  depth  close  to  these  critical  values. 

In  the  present  paper  we  construct  a  perturbation  solution  valid  at  the 
first  critical  value  of  the  depth.  Me  show  that  two  different  solutions  can 
exist  at  this  critical  value.  These  solutions  are  similar  to  the  "Milton 
ripples"  of  the  theory  of  gravity-capillary  progressive  waves  (Milton  (1915), 
Pierson  and  Fife  (1961),  Schwarts  and  Vandan-Broeck  (1979),  Chen  and  Saffaan 
(1979)). 

In  addition  we  use  the  numerical  scheme  derived  by  Vandan-Broeck  and 
Schwarts  (1981)  to  compute  the  solution  in  the  neighborhood  of  the  first 
critical  value  of  the  depth.  We  show  that  the  two  solutions  obtained  at  the 
critical  value  are  ambers  of  two  different  faailies  of  solutions. 

Me  formulate  the  problem  in  the  next  section.  The  main  results  obtained 
by  Concus  (1962)  are  summarised  in  Section  3.  The  perturbation  solution  valid 
at  the  first  critical  value  is  derived  in  Section  4.  The  numerical  results 
are  presented  in  Section  5. 


Sponsored  by  the  United  States  Army  under  Contract  Mo.  DAAG29-80-C-0041 .  This 
material  is  based  upon  work  supported  by  the  National  Science  Foundation  under 
Grant  Mo.  MC8-7927062,  Mod.  1". 


2.  romUXATXOM 


We  oouidir  the  t ln-ptr iodic  two-dinanalonal  potential  flew  of  a  fluid 
bounded  below  bp  e  horisontal  bottom  and  above  by  a  free  surface.  We  aeeuee 
the  notion  to  be  periodic  in  thf  horizontal  direction  with  wavelength  1.  We 
Measure  lengths  in  units  of  k”  ■  A/2*. 

Pol lowing  Coneus  ( 1962)  we  define  the  parameters  Y  and  6  bp  the 
relations 


Y 


(2.1) 


8 


1  ♦  Y 


(2.2? 


Here  o  is  the  surface  tenion.  For  4  «  1  the  capillary  effects  are  small, 
whereas  for  (1  -  4)  «  1  they  predominate . 

We  define  Cartesian  coordinates  such  that  the  action  is  syaawtrle  about 
the  vertical  plena  x  ■  0  and  such  that  y  ■  0  corresponds  to  the  naan 
level.  Let  k  h  denote  .the .mean  depth,  [kg(1  ♦  Y)]  '7  *  the  angular 
frequency,  [kg(1  +  Y)]  '  «  t  the  tins  and  a  the  amplitude  of  thf 

linearised  surface  wave  notion.  Then  we  define  *  *»  ak  and  let  Ck  n(x,t) 
denote  the  elevation  of  the  free  surface  above  the  scan  level  and 
c(g(1  +  Y)  1  V2k"3/>  the  velocity  potential. 

Zn  dimensionless  variables  the  action  of  the  fluid  is  described  by  the 
equations  (see  Coneus  (1962)) 


•*  0  in  0<x<*  and  -h  <  y  <  cn(x,t)  (2.3) 

(i  -  4)n  -  in^ii  +  "3/2  + 

♦  |  *<♦*  ♦  ♦*>  "  0  ««  Y  -  **»(x,t)  (2.4) 

♦y  -  w»)t  +  on  y  -  *n(x,t)  ,  (2.5) 

2f/0n  ■  0  on  x  -  0,  x  ■  *,  y  ■  -h  ,  (2.6) 

T)x  -  0  on  x  •  0,  x  ■  *  ,  (2.7) 


* 

/  *)(x,t)dx  -  0  , 
0 


(2.0) 


^♦(*»y»t  ♦  2*)  -  7#(x,y,t) 


(2.9) 


0*2* 

III  ♦(x,y,t)ain  t  eos  x  dt  dx  dy  -  0 
-h  0  0 
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(2.10) 


(2.11) 


v  ■  <*■  12  1/2 

•ad  III  ♦(*#y#t)  oo«  t  oos  x  dt  dx  dy  -  ~  *  (tanh  h)  '  • 

-b  0  0 

Aa  noted  by  Tadjbakhsh  and  Keller  (1960)  and  Oonous  (1962)  a  unique 
solution  does  not  exist  for  those  values  of  Yi-lat  which  tbs  frequency  of  the 
nth  spatial  harmonic  {nil  +  5(n  -  1 )tanh  nh)v*  is  an  integral  multiple  of 
the  fundamental  frequency  (tanh  h)*'2  this  yields  the  uniqueness  condition 

n[1  ♦  iiSL  -1I1W  I*  *  j2  for  {»  -  .  (2.12) 

tanh  h  J  ■  i**»«*. 

3.  PERTURBATION  SOLUTION  SATISFYING  THE  UNIQUENESS  CONDITION  (2.12) 

Following  Tadjbakhsh  and  Xaller  (I960),  Concur  (1962)  nought  a  solution 
as  an  expansion  in  powers  of  c.  Thus 


en  «  en°(x,t)  +  e2o1(x,t)  ♦  j  *3n2(x,t)  ♦  0(e4)  (3.1) 

*♦  -  e*°(x,y,t)  +  «V(x,y,t)  +  j  *392(x,y,t)  +  0(«4)2  (3.2) 

m  -  «0  ♦  et»1  +  i  *2m2  ♦  °(«3)  (3.3) 

The  solution  of  the  sero  >rder  solution  is  given  by 

'  °  ■  sin  t  oos  x  (3.4) 

9°  ■  (ttQ/sinh  h)cos  t  cos  x  cosh(y  +  h)  (3.5) 

w2  -  tanh  h  (3.6) 

This  solution  is  made  unique  by  imposing  the  condition  (2.12) 


Concus  (1962)  derived  the  following  expressions  for  the  first  and  second 
order  solutions 


2  -2 
Ml  ♦  m* 


_1  1  r  0  0  . 

n  -  e  Itht  + 


“2  -  -6 
M-  3M 

—  oos  2t]cos  2x 


1  -  36m! 


(3.7) 


♦1  ■  ®0  +  8  <W0  *  “03)t  “  16  (3“0  +  *l3)9in  21 

-  {3(m0  -  26m”3  -  (1  +  26)m”7]/16(1  -  36u”4)cosh  2h) 
sin  2t  oos  2x  oosh  2(y  ♦  h)  (3.8) 


StiAi.- 


«#1  -  0 

2 

n  -  b^ain  t  cos  x  +  b^3  sin  2t  cos  3x 
+  b^sin  3t  cos  x  +  b33«ln  3t  cos  3x 

m  &2  *  ®13°°*  t  008  3x  00811  3*y  + 

+  B^cos  3t  cos  x  co»h ( y+h )  +  B^cos  3t  oos  3x  cosh  3(y+h) 

.  -21*?  -  3(1+962)*»A  -  3(4+6d-9«2-2763)w"3  +  9{  H-56+4«2)n"7 

1  ^  0 _ 0 _ 0 _ 0_ 

32  (1  ♦  3«)(1  -  3«w“4) 


(3.9) 


(3.10) 


(3.11) 


]  (3.12) 


where  0  is  an  arbitrary  constant.  The  constants  bj^  and  B,  .  are 
defined  By  the  relations  (35)  and  (36)  given  by  Concus  (1962).  *3 

Tor  (  *  0  the  solution  (3.4)-(3.12)  reduces  to  the  solution  given  by 
Tadjbekhsh  and  Keller  (1960).  It  can  easily  be  checked  that  all  the  teres  are 
bounded  for  any  value  of  h  if  4  ■  0.  Thus  Tadjbakhsh  and  Keller's  solution 
is  a  satisfactory  third  order  solution  for  any  value  of  h. 

For  6  *  0,  sane  of  the  terns  appearing  in  u  ,  n  and  are 


r  ui  w  r  u  ,  ran  w  un  ww  mi  w.«  auiu 

unbounded  at  the  critical  values  of  depth  defined  by  the  relations 

1  -  3&IJ4  -  0 


(3.13) 


1  -  «(1  +  3n~)  -  0 


(3.14) 


These  critical  values  correspond  respectively  to  n  ■  j  -  2  and  n  »  j  «  3 
in  (2.12). 

In  the  next  eection  we  derive  a  perturbation  solution  valid  at  the  first 
critical  value  of  the  depth,  i.e.  at  the  value  of  the  depth  defined  by  (3.13). 


4.  PERTURBATION  SOLUTION  AT  THX  FIRST  CRITICAL  VALUB  OF  THX  DEPTH 

We  seek  a  perturbation  solution  of  the  fom  (3.1)-(3.3)  valid  when  (3.13) 
is  satisfied.  We  substitute  the  expansion  (3.1)-(3.3)  into  the  sys**n  of 
equations  (2.3)— ( 2.11)  and  collect  all  terns  of  like  powers  of  c.  The  terns 
with  c  to  the  first  power  in  (2.4)  and  (2.5)  are  given  by 


(1  -  «)n°  -  ♦  «*„♦£  -  0  on  y 

♦y  *  “0  on  y  «  0 


(4.1) 


(4.2) 


■nations  (2.3)  and  (2.6)-(2.11)  remain  unchanged  in  font  as  equations  for 
V  f  and  . 


Tha  t«ru  of  or dor  e2  in  (2.4),  (2.5)  and  (2.11)  oro  glass  by 

(1  -  5  )i)1  -  -  fQ  on  y  -  0 

(4.3) 

1 

♦J  “  -  «0  on  y  -  0 

(4.4) 

0  *  2*  - 

III  #1ooa  t  coo  x  dfc  dx  dy  ■  0 
-h  0  0 

(4.5) 

Boro  r0  and  G0  ora  dafinad  by 

F0  “  ‘  2  l(*x>2  +  ‘  V°*ty  "  “l*t 

(4.6) 

A  -  -°a°  -0.0  M  _0 

®o  ■  V,  -  n  ♦„  *  Vt 

(4.7) 

Aquations  (2.3)  and  (2.6)-(2.10)  ran* In  of  tha  Sana  fora  aa  aquations  in 
n,  41  and 

Tha  solution  of  tha  saro  ordar  problaa  dafinad  by  (2.3),  (4.1),  (4.2)  and 
(2*6)— (2-11)  is 

k 

i 

n°  “  sin  t  oos  x  +  A  cot  2t  oos  2x 

(4.8) 

•• 

*0  % 

•  ,'inh  hcoa  *  °°*  x  oosh(yfh)  -  2h*iP  2t  °°*  2x  °°*h  2(3T*b) 

(4.9) 

* 

h2  ■  tanh  h 

(4.10) 

Bara  A  is  aa  arbitrary  constant.  Thus  tha  solution  of  tha  saro  ordar 
solution  is  not  uniqua  whan  (3.13)  is  satisfiad.  . 

Diffarantiating  (4.3)  with  rsspset  to  t  and  substituting  n’  froa 
(4.4)  and  n V.  fron  (4.4),  aftar  diffarantiating  twloa  with  raapsot  to  x, 
wa  obtain 

«o «.  r-o 

(4.11) 

Bara  H0  is  dafinad  by 

Ho  "  V*  ♦  11  •  4,<s°  “ 

• 

(4.12) 

• 

8aparation  of  variables  yialds  for  tha  solution  of  (2.3) 
(2.6), 

subjact  to 

• 

♦'(X.y.t)  -  l  A  (t)oos  nx  cosh  n(y  ♦  h) 
n«o  n 

• 

(4.13) 

-5- 

Substituting  (4*13)  into  (4.11)  we  obtain 


whoosh  nh  A"(t)  +  ((1  -  «)n  ♦  6n3]sinh  nh  a  (t) 
u  n  n 


i  29 

■  —  f  H  cos  nx  dx 

V*  g  0 


(4.14) 


Bars  u  -  1  for  n  >  0  and  w  -  2  for  n  ■  0.  Using  (4.6)-(4.1Q)  vs  can 
rewrite  (4.14)  in  tha  form 

•JaJJU)  -  |  (3m3  +  »Q1)sin  2t  -  2A2 (m3 cotanh2 2h  +  3«3)sin  4t  (4.15) 

•Jooeh  h  A"(t)  +  sinh  hA^  [2*»1  ♦  J  (e"1  -  3«J  +  4«0ootanh  2h)]cos  t 


♦  |  [4*0cotanh  2h  ♦  e"1  +  21«3]coe  3t 


(4.16) 


«2oooh  2h  A£(t)  ♦  2(1  ♦  36) sinh  2h  A2(t) 

-  [»J  -  (1  +  2«)#»"1]  -  2Aa>1  -  6A»1«  -  4A«1  cotanh  2h}sin  2t  (4.17) 

«t2ooah  3h  A^(t)  +  (3  +  56) sinh  3h  Aj(t) 


AM.oos  t  72 

— -  ((4  +  486 ) cotanh  2b  -  (3  +  24«)t»"  -  3«*] 

AM  cos  3t 

- j - -  ((12  +  486) cotanh  2h  +  (3  ♦  246)n"2  -  21wJ) 


(4.18) 


«Qooah  4h  A” ( t)  +  (4  +  606)sinh  4h  A4(t) 

-  A2uQsin  4t((2  +  306) cotanh  2h  +  2«0cotanh22h  -  6mJ] 

e2oosh  nh  A"(t)  ♦  ((1  -  6)n  ♦  6n3]sinh  nh  A  (t)  -  0 
w  n  a 


(4.19) 


for  n  •  5,6,. 


(4.20) 


Prom  2.9)  and  (4.13)  it  follows  that  A_  must  ba  psriodic  in  t  with 
period  2*  for  n  >  1  and  from  (3.13)  and  (4.20)  that  A  -  0  for  a  >  5. 
The  periodicity  of  Aj  requires  the  coefficient  of  cos  tn  in  (4.16)  to  be 
equal  to  sero.  Thus 


1 


^  (3u3  -  a^  -  4utgCtft*nh  2h) 


(4.21) 


If  we  set  A  *  0  in  (4.15)-(4.21)  we  recover  the  system  of  equations 
derived  by  Concus  (1962)  for  the  first  order  solution.  In  particular  the 
solution  of  (4.17)  is  then  given  by 


3l»0  -  26wJJ3  -  (1  +  26)a“7] 
16(1  -  36a"4) cosh  2h 


sin  2t 


(4.22) 


This  solution  is  unbounded  since  (3.13)  is  assumed  to  be  satisfied.  Therefore 
we  do  not  set  A  -  0  in  (4.15). 

We  shall  determine  the  constant  A  in  such  a  way  that  the  solution  of 
(4.17)  is  bounded.  The  appropriate  compatability  condition  is  obtained  by 
multiplying  (4.17)  by  sin  2t,  integrating  with  respect  to  t  from  0  to 
2 *,  applying  integration  by  parts  twice  to  the  term  containing  A”(t)  and 
using  (3.13).  Thus  we  find  that  the  coefficient  of  sin  2t  in  the  right  hand 
side  of  (4.17)  must  be  equal  to  sero.  This  yields  the  relation 


3  [a?  -  (1  +  26)a"1] 

m  mi  ,  ■  U 

1  8  +  246  +  16UgCotanh  2h 


(4.22) 


Substituting  (4.21)  into  (4.22)  we  obtain 


,  3(«J  -  (1  +  26)a"1]  ' 

A  -  *  { - - - 2 - 5 - ^ - } 

[1  -  36  +  2aQcotanh  2h] [3»0  -  aQ  -  4aQcotanh  2h] 


(4.23) 


The  remaining  part  of  the  calculation  follows  closely  the  work  of 
Tadjbakhsh  and  Keller  (1960)  and  Concus  (1962).  Integrating  (4. 15)-(4. 19)  we 
obtain 

h0  “  “  16  (3"0  +  2t 


+  —•  (w0cotanh22h  +  3aQ)sin  4t  +  aflt  + 

At  -  -A[4«0cotanh  2h  +  a"1  +  21a3) [32sinh  h]"1 


A2  -  e^sin  2t 


(4.24) 


(4.25) 


(4.26) 


A3  "  Au^cos  t[(4  +  48*5)  cotanh  2h 

-  (3  +  246)<i>“2  -  3ti>2]  [( 12  +  206)sinh  3h  -  u2cosh  3h]-1 

-  AMqCos  3t [ ( 12  +  486 ) cotanh  2h 

+  <3  +  2*6)u~2-  21<i>2]  [ ( 12  +  20fi)sinh  3h  -  36«2cosh  3h]_1  (4.27 

A.  -  A2u  ain  4t((2  +  306) cotanh  2h  +  2o»Aootanh22h  -  6w2] 

4  0  0  0 

[(4  +  606)sinh  4h  -  16w2cosh  4h]“1  (4.28 

Hera  a  ,  0  and  a  are  constants  to  be  determined. 

Substituting  (4.13)  into  (4.3)  we  obtain 

1  1  4 

(1  -  6)n  -  6*1^  ”  PQ  “  “o  ^  A*  (t)cos  nx  cosh  nh  (4.29) 

n*0 

where  PQ  and  An(t)  are  defined  bp  (4.6)  and  (4.24)-(4.28).  The  function 
is  therefore  defined  as  the  solution  of  (4.29)  subject  to  (2.7). 

The  constant  a  in  (4.24)  is  evaluated  by  integrating  (4.29)  with 
respect  to  x  between  0  and  v  and  using  (2.7)  and  (2.8).  Thus  we  find 


1  a  a  - 

°0  "  8  “0  "  8  “o  +  ~2  (1  "  ootanh  2h>  (4.30) 

This  coatpletes  the  determination  of  the  first  order  solution.  2t  still 
contains  an  arbitrary  constant  o  .  This  constant  would  be  determined  at 
second  order  in  a  way  similar  to  the  way  A  was  determined  at  first  order. 
However  we  shall  not  do  this  in  this  paper. 

Equation  (4.23)  implies  the  existence  of  two  solutions  when  (3.13)  is 
satisfied.  Relations  (3.3)  and  (4.21)  show  that  one  solution  is  characterized 
by  a  frequency  larger  than  the  zero-order  frequency  and  the  other  solution  by 
a  frequency  smaller.  The  wave  profiles  given  by  these  two  possibilities  are 
illustrated  in  Figure  1.  These  solutions  are  very  similar  to  the  "Hilton 
ripples"  of  the  theory  of  gravity  capillary  progressive  waves  (Hilton  ( 1915) , 
Pierson  and  Fife  (1961),  Vanden-Broeck  and  Schwartz  (1979),  Chen  and  Saffman 
(1979)). 

In  the  next  section  we  show  that  these  two  solutions  are  members  of  two 
different  families  of  solution. 

5.  NUMERICAL  RESULTS 

Concus  (1962)  solution  is  satisfactory  for  values  of  the  depth  far  enough 
away  from  the  critical  values  (3.13)  and  (3.14).  The  solutions  derived  in 
Section  4  are  correct  at  the  critical  value  (3.13).  Perturbation  solutions 
▼slid  for  values  of  the  depth  near  but  not  equal  to  the  critical  value  (3.13) 


S  .  q 

•  «*>  *» 


Profiles  of  the  surface  of  the  standing  wave  at  t  “ 
curves  are  based  on  equation  (4.8)  with  6  ■  0.005  ai 
The  solid  curve  corresponds  to  A  >  0  and  the  broken 


could  be  obtained  by  using  the  P.L.K.  method.  An  example  of  such  a 
perturbation  calculation  can  be  found  in  Pierson  and  Fife  (1961). 

Xn  the  present  work/  we  compute  numerical  solutions  uniformly  valid  near 
the  first  critical  value  of  the  depth. 

Vanden-Broeck  and  Schwarts  (1981)  derived  a  numerical  scheme  4-o  compute 
pure  gravity  standing  waves.  Their  numerical  procedure  is  generalised  to 
include  the  effect  of  surface  tension  by  replacing  their  equation  (2)  by  the 
equation  (2.4).  The  numerical  procedure  then  follows  closely  the  method 
outlined  in  Section  XXX  of  their  paper. 

Numerical  values  of  m  as  a  function  of  36m”*  for  e  -  0.005  and 
h  ■  3  are  shown  in  Figure  2.  These  values  were  obtained  with  N  •  4  in  the 
equations  (15)  and  (16)  given  by  Vanden-Broeck  and  Schwarts  (1981).  Concus' 
perturbation  solution  for  m  is  represented  by  the  broken  line  in  Figure  2. 

It  is  defined  by  (3.3)/  (3.6)/  (3.9)  and  (3.12).  This  solution  is  unbounded 

when  36U-  -  1.  The  two  crosses  in  Figure  2  correspond  to  the  perturbation 

solution  of  Section  4.  They  are  defined  by  (3.3) ,  (3.6),  (4*22)  and  (4.23). 
These  two  solutions  are  in  fair  agreement  with  the  numerical  values. 

The  numerical  results  of  Figure  2  and  similar  results  obtained  for 
different  values  of  the  depth  indicate  that  the  solutions  derived  in  Section  4 
are  members  of  two  different  families  of  solutions.  One  family  of  solutions 
agrees  with  Concus  perturbation  solution  for  3 6m”  <  1  as  e  ♦  0  and  the 

other  family  agrees  with  Concus  perturbation  solution  for  36m”  >  1 •  Similar 

properties  were  found  by  Schwartz  and  Vanden-Broeck  (1979)  for  gravity¬ 
capillary  progressive  waves  in  the  neighborhood  of  the  first  critical  value  of 
the  capillary  number. 
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